Comparative motif discovery combined with comparative transcriptomics yields accurate targetome and enhancer predictions

This article is distributed exclusively by Cold Spring Harbor Laboratory Press for the first six months after the full-issue publication date. After six months, it is available under a Creative Commons License Attribution-NonCommercial 3.0 Unported License.-- et al.The identification of transcription factor binding sites, enhancers, and transcriptional target genes often relies on the integration of gene expression profiling and computational cis-regulatory sequence analysis. Methods for the prediction of cis-regulatory elements can take advantage of comparative genomics to increase signal-to-noise levels. However, gene expression data are usually derived from only one species. Here we investigate tissue-specific cross-species gene expression profiling by high-throughput sequencing, combined with cross-species motif discovery. First, we compared different methods for expression level quantification and cross-species integration using Tag-seq data. Using the optimal pipeline, we derived a set of genes with conserved expression during retinal determination across Drosophila melanogaster, Drosophila yakuba, and Drosophila virilis. These genes are enriched for binding sites of eye-related transcription factors including the zinc-finger Glass, a master regulator of photoreceptor differentiation. Validation of predicted Glass targets using RNA-seq in homozygous glass mutants confirms that the majority of our predictions are expressed downstream from Glass. Finally, we tested nine candidate enhancers by in vivo reporter assays and found eight of them to drive GFP in the eye disc, of which seven colocalize with the Glass protein, namely, scrt, chp, dpr10, CG6329, retn, Lim3, and dmrt99B. In conclusion, we show for the first time the combined use of cross-species expression profiling with cross-species motif discovery as a method to define a core developmental program, and we augment the candidate Glass targetome from a single known target gene, lozenge, to at least 62 conserved transcriptional targets.This work is funded by research grants from Research Foundation Flanders (FWO, grant G.0704.11N), University of Leuven (CREA/10/014 and PF/10/016), and Human Frontiers Science Program (RGY0070/2011). M.N.S. is funded by a PhD fellowship from FWO.Peer Reviewe

Neuronal cell types in the fly: single-cell anatomy meets single-cell genomics.

At around 150 000 neurons, the adult Drosophila melanogaster central nervous system is one of the largest species, for which a complete cellular catalogue is imminent. While numerically much simpler than mammalian brains, its complexity is still difficult to parse without grouping neurons into consistent types, which can number 1-1000 cells per hemisphere. We review how neuroanatomical and gene expression data are being used to discover neuronal types at scale. The correlation among multiple co-varying neuronal properties, including lineage, gene expression, morphology, connectivity, response properties and shared behavioral significance is essential to the definition of neuronal cell type. Initial studies comparing morphological and transcriptomic definitions of neuronal type suggest that these are highly consistent, but there is much to do to match these approaches brain-wide. Matched single-cell transcriptomic and morphological data provide an effective reference point to integrate other data types, including connectomics data. This will significantly enhance our ability to make functional predictions from brain wiring diagrams as well facilitating molecular genetic manipulation of neuronal types

dachshund Potentiates Hedgehog Signaling during Drosophila Retinogenesis

Proper organ patterning depends on a tight coordination between cell proliferation and differentiation. The patterning of Drosophila retina occurs both very fast and with high precision. This process is driven by the dynamic changes in signaling activity of the conserved Hedgehog (Hh) pathway, which coordinates cell fate determination, cell cycle and tissue morphogenesis. Here we show that during Drosophila retinogenesis, the retinal determination gene dachshund (dac) is not only a target of the Hh signaling pathway, but is also a modulator of its activity. Using developmental genetics techniques, we demonstrate that dac enhances Hh signaling by promoting the accumulation of the Gli transcription factor Cubitus interruptus (Ci) parallel to or downstream of fused. In the absence of dac, all Hh-mediated events associated to the morphogenetic furrow are delayed. One of the consequences is that, posterior to the furrow, dac- cells cannot activate a Roadkill-Cullin3 negative feedback loop that attenuates Hh signaling and which is necessary for retinal cells to continue normal differentiation. Therefore, dac is part of an essential positive feedback loop in the Hh pathway, guaranteeing the speed and the accuracy of Drosophila retinogenesis.MINECO Spain grants: (BFU2012-34324, BFU2015- 66040); Research Foundation—Flanders FWO grants: (G.0640.13, G.0791.14, PhD fellowship); Fundação para a Ciência e Tecnologia grant: (IF/01031/2012)

p53 shapes genome-wide and cell type-specific changes in microRNA expression during the human DNA damage response.

The human DNA damage response (DDR) triggers profound changes in gene expression, whose nature and regulation remain uncertain. Although certain micro-(mi)RNA species including miR34, miR-18, miR-16 and miR-143 have been implicated in the DDR, there is as yet no comprehensive description of genome-wide changes in the expression of miRNAs triggered by DNA breakage in human cells. We have used next-generation sequencing (NGS), combined with rigorous integrative computational analyses, to describe genome-wide changes in the expression of miRNAs during the human DDR. The changes affect 150 of 1523 miRNAs known in miRBase v18 from 4-24 h after the induction of DNA breakage, in cell-type dependent patterns. The regulatory regions of the most-highly regulated miRNA species are enriched in conserved binding sites for p53. Indeed, genome-wide changes in miRNA expression during the DDR are markedly altered in TP53-/- cells compared to otherwise isogenic controls. The expression levels of certain damage-induced, p53-regulated miRNAs in cancer samples correlate with patient survival. Our work reveals genome-wide and cell type-specific alterations in miRNA expression during the human DDR, which are regulated by the tumor suppressor protein p53. These findings provide a genomic resource to identify new molecules and mechanisms involved in the DDR, and to examine their role in tumor suppression and the clinical outcome of cancer patients

TOUCAN 2: the all-inclusive open source workbench for regulatory sequence analysis

We present the second and improved release of the TOUCAN workbench for cis-regulatory sequence analysis. TOUCAN implements and integrates fast state-of-the-art methods and strategies in gene regulation bioinformatics, including algorithms for comparative genomics and for the detection of cis-regulatory modules. This second release of TOUCAN has become open source and thereby carries the potential to evolve rapidly. The main goal of TOUCAN is to allow a user to come to testable hypotheses regarding the regulation of a gene or of a set of co-regulated genes. TOUCAN can be launched from this location:

Weakly Positioned Nucleosomes Enhance the Transcriptional Competency of Chromatin

Background: Transcription is affected by nucleosomal resistance against polymerase passage. In turn, nucleosomal resistance is determined by DNA sequence, histone chaperones and remodeling enzymes. The contributions of these factors are widely debated: one recent title claims ‘‘… DNA-encoded nucleosome organization…’’ while another title states that ‘‘histone-DNA interactions are not the major determinant of nucleosome positions.’’ These opposing conclusions were drawn from similar experiments analyzed by idealized methods. We attempt to resolve this controversy to reveal nucleosomal competency for transcription. Methodology/Principal Findings: To this end, we analyzed 26 in vivo, nonlinked, and in vitro genome-wide nucleosome maps/replicates by new, rigorous methods. Individual H2A nucleosomes are reconstituted inaccurately by transcription, chaperones and remodeling enzymes. At gene centers, weakly positioned nucleosome arrays facilitate rapid histone eviction and remodeling, easing polymerase passage. Fuzzy positioning is not due to artefacts. At the regional level, transcriptional competency is strongly influenced by intrinsic histone-DNA affinities. This is confirmed by reproducing the high in vivo occupancy of translated regions and the low occupancy of intergenic regions in reconstitutions from purified DNA and histones. Regional level occupancy patterns are protected from invading histones by nucleosome excluding sequences and barrier nucleosomes at gene boundaries and within genes. Conclusions/Significance: Dense arrays of weakly positioned nucleosomes appear to be necessary for transcription. Weak positioning at exons facilitates temporary remodeling, polymerase passage and hence the competency for transcription. At regional levels, the DNA sequence plays a major role in determining these features but positions of individual nucleosomes are typically modified by transcription, chaperones and enzymes. This competency is reduced at intergenic regions by sequence features, barrier nucleosomes, and proteins, preventing accessibility regulation of untargeted genes. This combination of DNA- and protein-influenced positioning regulates DNA accessibility and competence for regulatory protein binding and transcription. Interactive nucleosome displays are offered at http://chromatin.unl.edu/cgi-bin/skyline.cgi

Text-mining assisted regulatory annotation

Text-mining technologies can be integrated with genome annotation systems, increasing the availability of annotated cis-regulatory data

Cortical Neurogenesis Requires Bcl6-Mediated Transcriptional Repression of Multiple Self-Renewal-Promoting Extrinsic Pathways.

During neurogenesis, progenitors switch from self-renewal to differentiation through the interplay of intrinsic and extrinsic cues, but how these are integrated remains poorly understood. Here, we combine whole-genome transcriptional and epigenetic analyses with in vivo functional studies to demonstrate that Bcl6, a transcriptional repressor previously reported to promote cortical neurogenesis, acts as a driver of the neurogenic transition through direct silencing of a selective repertoire of genes belonging to multiple extrinsic pathways promoting self-renewal, most strikingly the Wnt pathway. At the molecular level, Bcl6 represses its targets through Sirt1 recruitment followed by histone deacetylation. Our data identify a molecular logic by which a single cell-intrinsic factor represses multiple extrinsic pathways that favor self-renewal, thereby ensuring robustness of neuronal fate transition